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Abstract 

The primary objective of this paper is to demonstrate 
the use of process-based manufacturing and assembly 
cost models in a traditional performance-focused 
multidisciplinary design and optimization process. 
The use of automated cost-performance analysis is an 
enabling technology that could bring realistic process- 
based manufacturing and assembly cost into 
multidisciplinary design and optimization. In this 
paper, we present a new methodology for 
incorporating process costing into a standard 
multidisciplinary design optimization process. 
Material, manufacturing processes, and assembly 
processes costs then could be used as the objective 
function for the optimization method. A case study 
involving forty-six different configurations of a 
simple wing is presented, indicating that a design 
based on performance criteria alone may not 
necessarily be the most affordable as far as 
manufacturing and assembly cost is concerned. 

Introduction 

The Multidisciplinary Design Optimization (MDO) 
methodology exploits the synergism of mutually 
interacting phenomena. The readers are referred to 
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recent review articles on MDO. \ 2 Traditional MDO 
tends to ignore cost and focuses primarily on vehicle 
performance criteria such as lift, drag, and range. If 
cost is included at all, then it is typically based solely 
on the weight of the vehicle. But this is inadequate 
and could even be misleading. High manufacturing 
cost could easily overwhelm any incentive to improve 
the design to the point of forcing the cancellation of 
the entire project. Determining the cost of 
manufacturing and assembly processes has been 
elusive in the past because of the difficulty of 
correctly modeling the cost of these processes. 

Typically the MDO processes focus on either 
optimizing the vehicle aerodynamic performance 3 or 
minimizing its structural weight. 4 5 The weight is 
indirectly related to the manufacturing cost, and the 
aerodynamic performance is related to operational 
cost. Both weight and performance play an important 
role in life-cycle cost. But they are not accurate for 
estimating the process-based manufacturing and 
assembly cost (PBMAC), which is directly related to 
the acquisition cost. Unfortunately it has been 
difficult to model the PBMAC in term of typical 
parameters and design variables used in a traditional 
MDO process. The purpose of this paper is to 
demonstrate the use of a PBMAC modeling tool with 
a performance analysis tool for cost-performance 
optimization. 

For our study, we have chosen to use the 
COSTRAN™ 6 code, 6 which is a commercial 
PBMAC. This code is an offshoot of a decade-long 
NASA effort 7 in developing PBMAC tools that is 


0The use of trademarks or names of manufactures in 
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traditionally used for aircraft trade study. The 
COSTRAN™ model is function of individual 
component parts such as spars, ribs, and skin, and it is 
a useful tool during the conceptual design phase of an 
aircraft. The goal of this paper is to demonstrate the 
use of commercial PBMAC in a traditional 
performance-focused MDO. The focus of this work is 
to determine the ’'what” (interface variables) and the 
"how" (interface methods) of integrating PBMAC 
tool with high-fidelity disciplinary models such as 
Finite Element Methods (FEM) structural models and 
Computational Fluid Dynamics (CFD) aerodynamics 
models. In the rest of this paper, the PBMAC model 
is first introduced. This will be followed by 
illustrative results obtained for the design of a generic 
wing. 



Assembly Cost Model (PBMAC 


systems with first-order velocity response to a step 
input as mathematically represented by the following 
equation: 

V = V 0 (l-e~%) (1) 

where V 0 is the steady-state process velocity, x the 
dynamic time constant, and t the process time. 

In general, t is governed by a major geometric 
property of the part, which could be its length, surface 
area, or volume. Using the terminology of reference 
15, this property is designated as X , the extensive 
variable for the process. 

The process velocity V can be equated to the first 
time derivative of X , i.e. V=dAAit. X can therefore be 
obtained by integration of V over time, resulting in 

A = V 0 [t-r(l-e^)] (2) 


The published literature abounds with 
articles and textbooks that advocate various PBMAC 
models. 811 A majority of these models rely on 
empirical data. In general, when manufacturing 
and/or assembly time is plotted against some design 
parameter on a log-log paper, a power law 
relationship between the variables can be determined. 
This procedure is the basis for a large number of Cost 
Estimating Relationships (CER) widely used in the 
industry. Another popular cost estimating procedure 
is the Response Surface Methodology (RSM) that 
relies primarily on multiple regression analysis. 12 
Finally the Genetic Algorithm (GA) is another cost 
estimating procedure tackling the problem from the 
standpoint of a biological phenomenon that enhances 
the successful processes while progressively 

eliminating the unsuccessful ones. 

All of the cost estimating methods 

mentioned above suffer from the following 
drawbacks: 1- Complete dependency on existing data, 
2- Application is limited to the range of available 
data, and 3- Unnecessary complication for early 

design optimization. Readers are referred to the 
literature for explanation of the drawbacks mentioned 
above. 

The work presented in this paper is 

supported by a commercial PBMAC. 6,7 The 

fundamental tenet of this PBMAC is a first order cost 
model first proposed in 1994. 15 This model was born 
out of an observation that many manual as well as 
automated processes can be represented as dynamic 


Equation 2 cannot be inverted explicitly for t. 
However two approximations can be made depending 
on the value of t relative to X such that: 


For t « x : t = J(2tX)/V 0 


b- 


For t » X : f = T + 



As suggested by Mabson (reported in reference 16), 
the above approximations can be combined into a 
single hyperbolic relation as followed: 


t = yla/V 0 ) 2 +(2tA/V 0 ) 


The validity of equation 3 can be seen in figure 1 
shown below. Other proofs are available in references 


14- 16. 


As indicated in reference 16, a total of 18 base time 
equations have been identified to directly relate the 
process time to the extensive variable under various 
conditions of operation. Bao provided a few case 
studies to illustrate the use of these equations. 17 

To illustrate the use of equation 3, consider the 
fabrication of a front spar for wing construction. 
Experience indicates that the V 0 and x values for a 
typical spar are respectively 2.4624 and 3.6934E+04. 
The extensive variable, X,, was determined to be the 
wetted area, i.e. area receiving machining, of the spar. 
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Therefore, if the spar’s wetted area is 100 in', then the 
fabrication time will be approximately 1732 minutes. 
Note that this fabrication time constitutes an overall 
time estimate without knowing all the details of part 
preparation, fabrication, and quality 
control/inspection requirements. During conceptual 
design phase, this time estimate is probably all that 
the designer needs to know for fabrication cost. 



Part Area, in2 


Figure 1- First-Order fit through industry estimates 
for abrasion operations (Reproduced from reference 
15). 


Preliminary Results 

For the purpose of demonstration, we have selected to 
use a generic wing, which is made of two spars, five 
ribs, and skin. Figure 2 shows the CAD representation 
of the generic model. The results are presented for 
two test cases: 1) cost comparisons for forty-six 
different concepts, and 2) cost optimization of generic 
wing concept. 



Figure 2 CAD representation of a generic wing. 


This model was parameterized using 
Multidisciplinary Aero/Structural Shape Optimization 
Using Deformation (MASSOUD 18 ) code. The 
MASSOUD code is based on a novel 
parameterization approach for complex shapes 
suitable for a multidisciplinary design optimization 
application. The approach consists of three basic 
concepts: 1) parameterizing the shape perturbations 
rather than the geometry itself, 2) utilizing Soft 
Object Animation (SOA) computer graphics 
algorithms, and 3) relating the deformation to 
aerodynamics shape design variables such as 
thickness, camber, twist, shear, and planform. 

The MASSOUD formulation is independent of grid 
topology, and that makes it suitable for a variety of 
analysis codes such as CFD and Computational 
Structural Mechanics (CSM). The analytical 
sensitivity derivatives are available for use in a 
gradient-based optimization. This algorithm is 
suitable for low-fidelity (e.g., linear aerodynamics 
and equivalent laminated plate structures) and high- 
fidelity analysis tools (e.g., nonlinear CFD and 
detailed Finite Element (FE) modeling). 


Figure 3 shows the parameterized model of a generic 
wing shown in Figure 2. This model has forty-five 
design variables, which consist of planform, twist, 
shear, camber, and thickness. 

Each set of forty-five design variables constitutes a 
design concept. All together, forty-six different 
design concepts were investigated. The basis for cost 
estimation per design concept is indicated in tables I 
and 2. 
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The interpretation of tables l and 2 should be as 
follows: the published values of V (l and T for an 
average spar were used as base values. The V 0 and x 
for all other wing components such as rib and skin 
were expressed in relative term compared to those of 
the base spar. Similarly the V 0 and x for the assembly 
of a typical wing were also used as base values. 
Values for the composite wing assembly were 
expressed in relative term compared to those of the 
aluminum wing assembly. It should be noted that 
wing assembly process should be separated from 
fabrication of skin, spar and rib because the former 
process depends critically on the perimeter while the 
latter process depends on the wetted area. Expressing 
all V 0 and x relative to those of the spar would be 
erroneous. Data in table 2 are representative of each 
of the indicated elements in a given year. 


Figure 3 Parameterized model of the generic wing. 


Table 1 : Basis for Cost Estimation of Generic 
Wing, V 9 and x 


Material: 

Aluminum 

V„ 

X 

Extensive 

Variable* 

Skin Fabrication 

1.228 

0.843 

A 

Rib Fabrication 

0.836 

1.122 

A 

Spar Fabrication 

1 

1 

A 

Wing Assembly 

1 

1 

B 

Material: 

Composite 

V 0 

X 

Extensive 

Variable 

Skin Fabrication 

0.871 

1.188 

A 

Rib Fabrication 

0.334 

0.280 

A 

Spar Fabrication 

0.588 

1.700 

A 

Wing Assembly 

0.714 

1.399 

B 


* Where, A is wetted area in inch 2 and 
B is perimeter in inch. 

Table 2: Basis for Cost Estimation of Generic 
Wing, Common Parameters 


Labor 

$60/Hour 

Material Cost: 


Skin 

$20/Lb 

- Rib 

$ 12/Lb 

Spar 

$ 15/Lb 

Fasteners 

$. 20/Unit 

Set Up and Delay Time per 

Not considered; 

operation 

Recurrence cost only 


For each design concept, the wetted areas for upper 
and lower skin, front and rear spar, and average rib 
were determined. Next, the perimeter for each of the 
above components was determined. Finally the data 
indicated in tables 1 and 2 were used to, first 
determine the fabrication cost of each component, 
second their assembly cost, and third and finally the 
total cost per design concept. Figure 4 shows the cost 
comparison for all forty-six different concepts, based 
on discrete choices of materials and shapes for a 
given structural topology, and given manufacturing 
and assembly processes. 

Figure 5 shows the cost comparison of individual cost 
factors for a given concept. 

For the first test case, i.e. aluminum wing, the 
parameterized model was embedded into an 
optimization process as shown in figure 6. 


Total Cost of Various Concept 
Aluminum Wing 


$212,000 

$211,000 

~ $ 210,000 
o 

w $ 209,000 
$ 208,000 


E 



Figure 4 Cost comparisons for forty-six different 
concepts. 
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Skin Fab - Skin Fab - 
Upper Lower 



Skins Assy Spars Assy 


□ Skin Fab - Upper 

■ Skin Fab - Lower 

□ Ribs Fab 

□ Spar Fab - Front 

■ Spar Fab - Rear 
B Spars Assy 

■ Skins Assy 

□ Ribs Assy 

■ Total Mtl Cost 


The optimization process is made of four modules: 
optimizer, geometry builder, cost estimator, and 
geometry constraints calculator. The optimization 
code CONMIN” was used for the optimizer module. 
As mentioned before, the MASSOUD code was used 
to parameterize the geometry. The cost estimating 
concept described previously was used to estimate the 
cost of a generic wing. The total wetted skin, rib, and 
spar areas were constrained to stay below the baseline 
design. 

Figure 7 shows preliminary optimization result for 
the generic wing shown in figure 2. The cost was 
reduced by more than 1.8%. 


Figure 5. Cost comparisons of individual cost factors 
for a generic wing. 



Figure 6. Optimization Process 



Figure 7. Cost optimization. 


Discussion 


Cost consideration is among the most important 
elements in any multi-disciplinary design 
optimization scheme. There are many kinds of cost 
involved in a typical airplane program. As described 
by Roskam, 20 there are costs associated with the 
planning and conceptual design, with preliminary 
design and system integration, with detail design and 
development, with manufacturing and acquisition. 
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design and system integration, with detail design and 
development, with manufacturing and acquisition, 
with operation and support, and with disposal. This 
paper deals strictly with the first type of costs, notably 
costs associated with the planning and conceptual 
design. As indicated earlier, the MDO community so 
far tends to treat cost as solely based on the weight of 
the vehicle. The case studies included in this paper 
indicate that fabrication and assembly costs are much 
more significant than material costs - as expressed by 
weight- and should be part of the optimization 
scheme. 

Even at the conceptual design phase, there is a need 
to incorporate the costs of fabrication and assembly 
of the major components such as spars, ribs, and 
skins. Using the first design configuration as a typical 
design, the following table reveals how dominating 
fabrication and assembly costs were over material 
costs. 



Mtl 

Mfg 

Assy 

Total 

Wing 

Front 

Spar 

5.5% 

21.3% 

73.2% 

12.3% 

Rear 

Spar 

4.3% 

19.5% 

76.2% 

11.4% 

5 Ribs 

3.9% 

23.4% 

72.7% 

39.9% 

Upper 

Skin 

7.5% 

33.4% 

59.1% 

18.4% 

Lower 

Skin 

7.6% 

33.9% 

58.5% 

18.0% 

Total 

Wing 

5.5% 

26.4% 

68.1% 

100% 


From the above percentage table, it can be said that, 
in general material cost was only about 5% of the cost 
of fabrication and assembly. Also, fabrication cost of 
either spar or rib was about 30% of corresponding 
assembly cost, while fabrication cost of skin was 
about 50% of assembly cost. The numbers quoted 
above are close to industry standards. 

As to the cost comparison of the forty-six different 
design concepts, while the magnitude of the overall 
cost reduction was less than 2%, the point was that 
the proposed cost model was detailed enough to 
accommodate all design concepts. Furthermore it 
could be easily incorporated in any multi-disciplinary 
optimization methodology. 


Conclusions 


We have demonstrated the use of process- based 
manufacturing and assembly cost models in a 
traditional performance-focused multidisciplinary 
design and optimization process. Three major 
conclusions can be drawn from this paper. First the 
weight may not be directly related to cost, and 
minimizing the weight without considering the 
manufacturing and assembly costs may increase the 
overall cost. Second the analytical cost models can be 
incorporated in a traditional MDO process. And third, 
the fabrication and assembly costs could drive the 
optimization process to minimize the actual cost of 
the part being considered. 
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